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Abstract 



> 

o 

OO ' We prove the existence of a class of plane-symmetric perfect fluid 

. cosmologies with a Kasner-like singularity of cigar type. These so- 

lutions of the Einstein equations depend on two smooth functions of 
one space coordinate. They are constructed by solving a symmetric 
hyperbolic system of Fuchsian equations. 



99' 1 Introduction 

^ , The Bianchi I spacetimes are the simplest class of anisotropic cosmological 

^ ' models. They admit a 3-parameter abelian group of isometries acting on 



spacelike hypersurfaces [17|. The dynamics of solutions to Einstein's equa- 
tions with this symmetry type are particularly well-understood in the case of 
a perfect fluid as matter source. For a physically reasonable, non-stiff equa- 
tion of state and with an appropriate choice of time-orientation all solutions 
isotropise in the future and in the past any solution which is not FRW must 
become Kasner-like as the initial singularity is approached. Solutions with 
all values of the Kasner exponents are realised [^] . 

Under the additional assumptions of local rotational symmetry and re- 
flection symmetry much can also be said about the asymptotic behaviour 
of collisionless gas cosmologies (i.e. solutions of the Einstein- Vlasov equa- 
tions). Once again all solutions isotropise in the future. For any solution 
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which is not FRW the asymptotic generaUsed Kasner exponents (i.e. the 
eigenvalues of the second fundamental form divided by the mean curvature) 
in the past are either (—5, |, |) or (0, ^, ^) ||l5|. The first case is Kasner like 
and is generic in that it occurs for an open dense set of initial data. Note 
that no solution of the (0, ^, ^) type occurs in the fluid case and that the 
Kasner exponents (1,0,0) are not realised with the collisionless gas. 

It is natural to ask whether the above qualitative features are main- 
tained if one relaxes the requirement of spatial homogeneity somewhat. 
Certainly there are many inhomogeneous FRW-like solutions of both the 
Einstein-perfect fluid equations and the Einstein- Vlasov equations, as has 



been shown in the framework of isotropic singularites |l2|]-[^, y], Kas 



ner like solutions of the Einstein- Vlasov equations have been investigated 



by Rein [12| under the assumption of plane symmetry, which is the simplest 
inhomogeneous generalisation of LRS Bianchi I symmetry. By proving a 
global existence theorem he is able to show that an open set of C°° data 
on a regular Cauchy surface evolves back to a Kasner like singularity of 
the (— 1,|,|) type, agreeing with the picture described in In the 

present paper we study plane-symmetric, non-stiff 7-law perfect fluid cos- 
mologies with asymptotic geometry of this same type. The results obtained 
can be thought of as providing some hard evidence for certain aspects of the 
general picture of cosmological singularities proposed by Belinskii, Khalat- 
nikov and Lifschitz (BKL) many years ago |^. This picture of the generic 
solution of Einstein's equations, based on formal calculations, has a number 
of key elements: 

(i) Near a singularity the evolution at different spatial points should de- 
couple so that spacetime is approximated by a spatially homogeneous 
model at each point. 

(ii) The matter content of spacetime should in general not have a signifi- 
cant effect on the dynamics near the singularity. Thus the dynamics 
should be well approximated by vacuum spacetimes. 

(iii) The most general spatially homogeneous cosmology is modelled by 
a spacetime of Bianchi type VIII or IX and vacuum spacetimes of 
Bianchi types VIII and IX typically exhibit Mixmaster behaviour. 

One should note that there are exceptional models which do not fit the 
BKL picture. Spacetimes with isotropic singularities are an example, but 
they depend on fewer free functions than the general solution and so do not 
constitute a counterexample. Also, while the BKL picture predicts in general 
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complicated, oscillatory Mixmaster behaviour one may still investigate (i) 
and (ii) for less complicated Bianchi type I-like models. Rigorous results 
for Gowdy spacetimes have been obtained by Isenberg and Moncrief ||^] 
and Kichenassamy and Rendall Recall that the Gowdy spacetimes are 
defined as spatially compact solutions of the vacuum Einstein equations 
which have two commuting Killing vectors and in which the so-called 'twist 
constants' vanish. In the so-called Fuchsian method is used to construct a 
class of Gowdy spacetimes which have Kasner-like singularities and depend 
on the maximum number of free- functions. In the authors specialize to so- 
called polarized Gowdy symmetry and use energy estimates to show that all 
solutions are 'Asymptotically Velocity Term Dominated' near a cosmological 
singularity. This is one way of saying that the field equations spatially 
decouple as in (i). 

The method used in this paper to investigate perfect-fluid spacetimes is 
very similar to that employed in It consists in taking a two-parameter 
family of exact Bianchi I solutions, making an inhomogeneous perturba- 
tion and solving a singular (Fuchsian) system of field equations for the per- 
turbation near the singularity 't = 0'. In this way we construct a two 
free-function family of C°° inhomogeneous cosmologies with the required 
asymptotic Bianchi I like behaviour. This family is 'generic' in the sense 
that the general plane-symmetric perfect-fluid spacetime depends on just 
two free- functions. 

The main result is stated formally as Theorem 6.1 at the end of section 

6. 

In section 7 we show that the spacetimes constructed have crushing sin- 
gularities. 

2 Exact solutions 

The general Bianchi I solution of the Einstein-perfect fluid equations, as 



given in |16| is 

ds^ = -B'^-'-'^dT^ + T^P^B'^^dx' + T^P^B'^^dy' + T^P^B^'i^dz^ (1) 
where 



5^-^ = a + m^T'^-^ a > 0, m > (2 

Pi + P2 + P3 = 1 , Pi + P2 + Ps = 1 ' qi = '^-Pi (3 
(The case a = is an FRW model and is excluded in what follows.) 
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The density of the fluid is given by 

and 7 is the polytropic index, so that the pressure is given by P = (7 — 
It is a consequence of the momentum constraint that in Bianchi type I 

the fluid cannot tilt and hence the 4- velocity points in the direction of 
Now if one wants to consider the plane-symmetric cases of (1) then there 

are two choices for the Kasner exponents pj, namely (— ^, |, |) and (1, 0, 0). 

In this paper we will restrict attention to the (— |, |, |) case and consider 

inhomogeneous perturbations of the metric which depend on a coordinate 

in the '— | direction'. 

When we come to write down Einstein's equations for a plane symmetric 

metric we will use conformal coordinates, i.e. coordinates for which the 

metric takes the form 

ds^ = e2^(*'^)(-di2 + dx^) + ji(^t, x){dy^ + dz^) (5) 

The model solution (1) in the case pi = (— |, |, |) may be written in 
conformal coordinates by making the transformation 



T ^ t = / 7^ — ds (6) 

Then the metric (1) takes the form (5) with 

A = -^logT + -^log(a + mV-T) , R = T^/^ (7) 
o ^ ~ 7 

and r given implicitly as a function of t by the relation (6) . 

3 The Einstein-perfect fluid equations 

We will assume that spacetime is filled with a polytropic perfect fiuid, so 
that the stress tensor takes the form 

T"^ = (p + P)u''u^ + Pg""^ (8) 

with u"ua = -1 and P = (7 - l)p, 1 < 7 < 2. 

With this stress tensor the Einstein evolution equations for the metric 
(5) are 

Rtt - Rxx = (2 - j)Re^^p (9) 
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Att - - \r-\RI - Rl) = -\ie'^P (10) 
while the constraints read 

R-^Rxx-\R~''{Rl+Rl)-R~^RtAt-R-^RxAx = -e''^p{jiv^f + l) (11) 

- R-^Rtx + R-\RtAx + AtRx) + ^R-^RtRx = -e''^^pv%^ (12) 

where = e'^u^ and hence (v^)"^ — (f^)^ = 1. 

The Euler equations VTq,^ = take the following explicit form 

(1 + ^{v^f)pt + 2-fp{v^)t - 2-ipAt{l + {v^f) + ^v%^px 

+ 7P {At{?>{vy + {v^f) + v''v\AAx + R-'Rx) + [v^fR-^Rt] = (13) 
and 

^v^v'pt + ^p I (^ (^^)' + (^°)' ^ ((,1)^ _ A,v^^ _ ^^^1 

+ (7(1 + {v^f) -l)px + 2-ipv\v^)x 

+ 7p {v^v\2At + i?-^i?t) + (v^fR-^Rx + (1 + 2(i;^)2)yl^ + 2^;%^^^} = 

(14) 

4 Inhomogeneous perturbations 

We now choose two strictly positive smooth functions a{x), m{x) and seek 
solutions p, v^, A, R of (9)-(14) which approach the model forms (4)- 
(7) as t ^ 0. 

To be specific we make the following ansatz 

^ = -^logr + --^log(a(x) +m^(x)r^~'^) +tA (15) 
3 2 — 7 

R = T^I^{l+R) (16) 
logp = log/x + tt(T-|)+^0 (17) 
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where < e < |(2 — 7), r(t, x) is given by (6) and 



The aim is to find solutions of the Einstein-perfect fluid equations for 
which 

A, R, (j) and tp tend to zero as t tends to zero. The choice of V"*^ ensures 
that the momentum constraint is satisfied at t = (see section 6). 

It is convenient to write the field equations for A, R, (() and il) in first 
order form. To do this we introduce the following new variables: 

U = Rt, Q = dt{tA), X = R^, Y = tA^, S = t-^R 

In terms of these variables the evolution equations (9)- (10) take the form 

tXt = tU^ (19) 

tYt = tQ^ (20) 

tRt = tU (21) 

tAt + A-Q = (22) 

tdtS + S-U = (23) 

tUt + 2U = tX^ + (^-2ti^^-^^F{t,x) - ^mhrl-^^ U 

+ + R){2- 7)(a + mV2-^)(e2*^+*'+^""'^ - 1) (24) 

tQ, - ic/ = tY, + tAl, - ^jmHr-l-^a + m\^-^){e^'^+''^'-''^ - 1) 
+ lt(l + Ry (u^ + ^ri-^m^l + R)U + ((logr4/=^)^(l + R) + Xf^ 

+ |tl(2-T)F(t, x) + i(l + p(2-7)^(i^ (25) 

where AP stands for the first two terms in (15), the function F{t,x), contin- 
uous in t and smooth in x is given by the relation 

lair-t = 1 + t^^^-'^^F{t, x) (26) 
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The Euler equations may, after some rearrangement, be written 

+2i2-2f^^, + (1 + 2t2-2fV;2)(^0 + y) + (logr4/3), 
+(1 + R)-^X + ii-PV^^0(4ri-^ + 4Q + (1 + i2)"^C/} 

-7 |q + ^t'^-'^P^^^ar-^/^ + 4mV2-^ + 3Q) + Av^t^-^iljiAl + F) 

^ (27) 



and 



-t^-W |4mVi-^ + + 4Q + (logr^/^)^ ^ 

-72t2-2p^2^-|(^^^2^2-7^+^^i-P^(2ti-J'Vt+(2-2p)t-PV+t'°((log/x)^+i^+^-^'0^)) 
+7 (q + 2i2-2fV'(-^«r-t + ^mVi-^ + Q) + 47;0ti-PV(^S + 

+^ (^(^1±!!^^ [/ + tr-W-^v^^a + mV-^) + ^^^t^-f^ ((logri), + J^)) } 

+5 / 2(1 + w^-^t-^'v + i^(i + t'-'^v')^ + — ^ 

- 27V^t2~2p^^ _ ^ 2i2-2p^2)(^o ^ Y) (28) 
where p = |(2 — 7) and = + t^'i/i. 



7 



5 Existence and uniqueness of solutions 



The Einstein evolution equations (19)-(25) take the form 

tai{wi)t + tbi{wi)x + ciwi =t^Gi{wi,(f),t,x) (29) 

where wi stands for U, A, Q, X, Y,R,S, ai is the 6x6 identity matrix, 6i is 
a constant symmetric matrix and ci is a constant matrix with eigenval- 
ues 0, 1, 2. Gi is a vector- valued function, continuous in t and smooth in all 
other arguments. 6 is a strictly positive constant. 

The Euler equations (27)- (28) may be symmetrised by multiplying one 
of the equations through by the smooth function of t^~P'tp given by 

Vi + 2t2-2p^2 yo yv 7 i + 7t2-2p^2 J 

(30) 

In this way (27)- (28) give rise to a system of the form 

C2W2 = t^G2 {wi ,W2,t,x) (31) 

where W2 stands for cp, ^, 02 is a diagonal positive definite matrix, smooth 
as a function of t^~Pilj (and equal to the identity at t^~'Pil) = 0), 62 is a 
smooth symmetric matrix and C2 is diagonal with positive eigenvalues. G2 is 
a vector- valued function, continuous in t and smooth in its other arguments. 
Combining (29)-(31) we see that the field equations may be expressed as 

tai{t^-Pilj){wz)t + tbsii/j, t){w3)^ + C3W3 = t^Gs {W3, t, X, ) (32) 

where 03 is smooth and positive definite, 63 is smooth and symmetric and C3 is 
constant with positive eigenvalues. G3 is a vector-valued function, continu- 
ous in t and smooth in all other arguments. 

Now note that we may choose ^ = ^ and | (7 — 2/3) = q/n for some 
positive integer n and some q > I wlog. With this choice of 5 it is then 
convenient to make the change of time coordinate t ^ s = tn . Equation 
(32) then takes the form 

a3{w3)s + b3{w3)x + ^C3W3 = Gs (33) 

where as = n~^a3 and 63 = s"~^63. 

It is also convenient to multiply the right hand side of (33) by a smooth 
cut-off function of x with compact support. Call the result G3. 
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We are aiming to solve the equation 

03(^3)5 + h(.W3)a: + -C3'U;3 = G3 (34) 

s 

with ^3(0) = on [0,T] x M. A solution of this equation will be a local (in 
space) solution of (33) by the finite speed of propagation. 

For the application of energy estimates to a Fuchsian system it is impor- 
tant that the singular part of the system be positive definite. Unfortunately 
the matrix C3 in (34) is not positive definite even though its eigenvalues are 
all positive. The situation can however be remedied by making a change of 
variables. First, because C3 has positive eigenvalues one can, by Lemma 2.1 
of [^8| obtain an approximate solution Ui of (34) which satifies this equation 
up to 0(s*) for any positive integer i, i.e. 

sa3{ui)s + sb3{ui)^ + C3Ui - sGs = 0{s') 

as s — > 0. The function Ui{s,x) is continuous in s at s = and smooth 
for s > 0, with nj(0, x) = 0. 

Now define the new variable W3 according to W3 = s~^{w3 — Ui). Then 
as an equation for (34) becomes 

a3{w3)s + b3{w3)x + -C3W3 = G3 

s 

where £3 = C3 + il and G3 has the same properties as G3. 

We now assume that i was chosen large enough to make 03 positive defi- 
nite and in the following we drop the tildes from W3 and C3 for convenience. 

Now G3 contains terms like fit, x)g{w3)w3 or just h{t, x) where / and h are 
continuous in t and smooth with compact support in x and g is smooth. We 
approximate such functions /, h uniformly in the Sobolev space H'^{W) on [0, T] by 
sequences of smooth (C°°) functions of compact support /i" (see Lemma 
A.l of the Appendix and see Q for an introduction to Sobolev spaces). 

Consider the equation 

~a3{s.W-3){^-,+% + b3{s,W^)iw^+\ + ^C3<+1 = E/^S^KX"'' + K 

(35) 

where the right hand side represents an approximation to G3. Since this 
equation is linear, all coefficients are smooth and C3 is positive definite, we 
can iteratively obtain smooth solutions ^3 on [0, T] with 1^3(0) = (0) by 
Theorem 9.1 of Lemma 5.1 The sequence is uniformly bounded 
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in H'P and for the element ip"^ of w"^ we have that sdsip^ is uniformly bounded 
in HP^^ on [0, Ti] for some Ti <T. { p is chosen large enough to apply the 
Sobolev embedding lemma wherever needed). 

Proof We use induction. 
Suppose 

(i) WwsWpis) <Riors<T 

(ii) ||s9^V'"||p-i <LioT:s<T 

Apply X derivatives of order a to (35) to get 

ds{w^)^sV''w'^^+^ + fe3+^a^V"u'^+^ + Jc3V°u;^+^ = M° (36) 

where contains terms like 

asV" [d^%d,w^+') (37) 

and the terms coming from the /" and /i". 

Now take the inner product of (36) with V^Wg"*"^ and use the symmetry 
of 63, the boundedness of sdgip"' and the uniform equivalence of the norm 
of V"w?+^ to 



a3V"u;^+^ • V°u;^+^ dx (38) 
in a standard way |11], to get the following inequality 

\\p{s)<Kil \\w'^+^\\p + \\wl\\p + K2dr (39) 

JO 



,,n+l I 



(We used the fact that the /" are bounded in and we also used the posi- 
tive definiteness of C3 to throw away the singular ^ terms which occur in (36)) 
If Ti is now chosen sufficiently small Gronwall's inequality gives Huig llp(s) ^ 
for s < Ti. 

Prom the differential equation (35) one may easily use a Moser estimate 
[11| to obtain ||s(9sV"''''^llp-i ^ L for s <Ti ii L was chosen large enough. 



Hence the lemma. 

Note: From the form of 03 the bound on sdsip^ is enough to give control 
over the time derivative of (13 (wg) which occurs in the right hand side of the 
energy estimate for w^^^ . 

Lemma 5.2 The sequence w"^ converges uniformly in on [0, T2] for 
some T2 <Ti. 

Proof. By subtracting the equation satisfied by from that satisfied 

"3-' 



by Wo ^ and again exploiting the positive definiteness of C3 one obtains the 
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following energy estimate 



- W3~%i^) <K r {1^3 - + Ik""' - «^""'ll2 

^0 



+ iir-'-r-l2} (40) 

(Note; we need boundedness of the quantity (03(^3) — d3{w^~^))dsW^ to 
obtain this inequality. But this follows from the form of 03 and Lemma 5.1) 
Put a" = - w^-^h, 7" = 11/""^ - r~^\\2- Then we have 



Gronwall implies 
which implies 

Thus 



a"(s) <K r [a"" + a"-i + 7"] dr (41) 
Jo 

a"" <K r [a"-i + 7"] dr (42) 
Jo 

sup a"" <K [ sup a"~^ + 7" (43) 

s<T2 Jo S<T2 



sup a" < KT2 sup a"-^ + sup 7" (44) 

s<T2 s<T2 s<T2 

We assume that the sequence /„ was chosen to converge so fast that X] 7" < 
00. We also choose T2 so that KT2 < 1. It now follows by the argument of 
Corollary 5.4 in that tUg is uniformly Cauchy in on [0,T2]. Hence 
the lemma. 

The inequality of Gagliardo-Nirenberg [11| now shows that w'^ converges 
uniformly in HP~^ to some w on [0,r2]. p was arbitrarily large. From 
the differential equation one see that quantities dgw'^ converge uniformly 
on (e, T2] for e > 0. Hence w is by the Sobolev imbedding theorem. 
A standard continuation argument |14| for regular symmetric hyperbolic 
systems can be applied away from s = to get that w is C°° for s > 0. 

Suppose now that we have two solutions w, w of the field equations, 
smooth for s > 0, continuous at s = and such that w{0) = w{0). Then by 
taking the difference of the equations satisfied by w and w one may obtain 
the following energy estimate. 

ds\\w — w\\2 < K\\w — w\\2 (45) 

for s > 0. (Again one needs control over {d^{w) — d3{w))dsW, which follows 
from the equation) 



11 



Hence 

\\w — w\\2{s) — Ww — li^lll(e) < K J \\w — w\\2 dt (46) 

which impUes 

\\w - w\\l{s) < K f \\w - w\\l dt + \\w - w\\l{€) (47) 
Jo 

Gronwall's inequality now implies 

\\w -w\\l{s) <e^'\\w -w\\l{e) (48) 
Since e > was arbitrary, we must have w = w. 



6 The constraints 

Define constraint quantities Co, Ci by 

Co = R-'Rxx - + Rl) - R-^RtAt - R-^R^A^ + (?^p{\ + liv^ f) 

(49) 

Ci = -R-^Rtx + R-^iRtA,, + ^ti?^) + ^R-^RtR^ + e^^jpv%^ (50) 

If the evolution equations (9)-(10), (13)-(14) are satisfied then a calculation 
shows that the following hold 

dtCo = -d,Ci - ^Co - ^Ci (51) 

dtCi = -d,Co - ^Ci - ^Co (52) 

One also calculates that the quantities Cq = tCo, Ci = tC\ tend to zero 
as t tends to zero. These quantities satisfy the following 

tdtCo + - l) Co = -td^Ci - (53) 

tdtC^ + - l) Ci = -td.Co - ^Co (54) 

Now R:c/R is 0(1) and tRt/R = 1 + 0{t^) for some 5 > 0. It thus follows 
from the ideas of section (5) that Co and C\ are identically zero and hence 
the constraints are satisfied. 
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The results of sections 5 and 6 can now be summarised as follows: 



Theorem 6.1. Given two smooth, strictly positive functions a{x), m{x) on M 
there exists a unique solution {g^^y^p^u'^) of the Einstein equations coupled 
to a 7-law perfect fluid on x (0, T) satisfying 

ds^ = e2^(*'^) {-df + dx^) + R{t, x){dy^ + dz^) 

and 

A = A^ + o{t) 
i? = r5(l + o(l)) 

log/9 = log// + o(t^-i(2-^)+^) , 0<e<-(2-7) 

tdtA = tdtA^ + o{l) = ~+o{l) 
dtR =^a + o{t) 

where 1 < 7 < 2 and 

A^ = -^logT + -^log(a(x) + m2(x)r2-^) (55) 
3 2 — 7 

. r «^ ^ 34 

Jq a + m'^s'^ ' 4a 

7 Crushing singularities 

The extrinsic curvature of the hypersurfaces of constant t is given by 

K = e-^\-\a + m^T^-^)^{a + 2m^T^-^ + t^/^{Q + U{1 + i?)"^)) (56) 

It follows that in the limit as t ^ we have K ~ a^'y+'^y^'^-'^h'i . Thus 
if a is chosen uniformly positive we have a crushing singularity at t = 0. 
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8 Conclusion 



We have proved the existence of a family of plane-symmetric perfect-fluid 
cosmologies which have a (— |, |, |) Kasner-like singularity and which de- 
pend on the full number of free functions for perfect-fluid spacetimes with 
this symmetry type, namely two. This supports the idea of Belinskii, Kha- 
latnikov and Lifschitz that the Einstein equations should generically decou- 
ple spatially near a cosmological singularity and that the matter content of 
spacetime should have a negligible effect on the dynamics near the singular- 
ity. 

It would be interesting to know if the result obtained in this paper can 
be extended to perfect-fluid spacetimes with polarized Gowdy symmetry (i. 
e. diagonal G2 cosmologies). The formal expansions of BKL and the work 
of Isenberg and Moncrief suggest that this should indeed be the case. It 
may well be that the Puchsian algorithm can once again be applied in this 
more general setting, but this has not yet been attempted. 



A Approximation of continuous functions by smooth 
functions 



Lemma A.l Suppose we are given a function F{t, x) on [0, T] x R, which 
is continuous in t and smooth in x. Suppose also that F has fixed compact 
support in x. Then there is a sequence of smooth functions F^it, x) on [0, T] X 
M such that F„ ^ F in W{R), uniformly on [0,r]. 

Proof First define 

' F{t,x) 0<t<T 
F{t,x) = }. F{0,x) -l<t<0 (57) 
F{T,x) t>T 



and let 



Fnit,x)= / Ji{t- s)F{s,x) ds 



(58) 



where is a smooth function, supported on (— e, e), such that J Je{x)dx = 1. 
Then 



\F{t,x)-Fnit,x) 



Ji (i - s)iF{t,x) - F{s,x)) ds 



< sup \F{t,x) - F{s,x)\ 
l*-s|<- 
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Now 



sup / \{F - Fn){t,x)\^ dx < / sup ( sup \F{t,x) - F{s,x)\) dx 

0<t<TjR JRO<t<T y\t-s\<^ J 

(59) 

Clearly the integrand in the right hand side of (59) tends to zero for each x by 
the uniform continuity of F in t. It is also clear that the sequence of inte- 
grands is bounded by an integrable function. Thus Fn F in uniformly 
on [0, r] by Lebesgue's dominated convergence theorem. One may now ap- 
ply the above argument to partial derivatives to get the desired result. 
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